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The harmonic oscillator is one of the simplest physical 
systems but also one of the most fundamental. It is ubiq- 
uitous in nature, often serving as an approximation for a 
more complicated system or as a building block in larger 
models. Realizations of harmonic oscillators in the quan- 
tum regime include electromagnetic fields in a cavity 0]- 
O and the mechanical modes of a trapped atom [4j or 
macroscopic solid |5 |. Quantized interaction between two 
motional modes of an individual trapped ion has been 
achieved by coupling through optical fields |6|, and en- 
tangled motion of two ions in separate locations has been 
accomplished indirectly through their internal states (7). 
However, direct controllable coupling between quantized 
mechanical oscillators held in separate locations has not 
been realized previously. Here we implement such cou- 
pling through the mutual Coulomb interaction of two ions 
held in trapping potentials separated by 40 /im (similar 
work is reported in a related paper (8)). By tuning the 
confining wells into resonance, energy is exchanged be- 
tween the ions at the quantum level, establishing that di- 
rect coherent motional coupling is possible for separately 
trapped ions. The system demonstrates a building block 
for quantum information processing and quantum simu- 
lation. More broadly, this work is a natural precursor to 
experiments in hybrid quantum systems, such as coupling 
a trapped ion to a quantized macroscopic mechanical or 
electrical oscillator I9HT3I. 

The direct coupling of atomic ions in separate potential 
wells is a key feature of proposals to implement quantum 
simulation fT4ifT6lL and it could allow logic operations to 
be performed in a multi-zone quantum information proces- 
sor (lOl [TTl [18) without the requirement of bringing the ion 
qubits into the same trapping potential. Moreover, the cou- 
pling could prove useful for metrology and sensing. For ex- 
ample, it could extend the capabilities of quantum-logic spec- 
troscopy @ Q3D [20| to ions that cannot be trapped within 
the same potential well as the measurement ion, such as op- 
positely charged ions or even anti-matter particles (9l [Toll . 
Coupling could be obtained either through mutually shared 
electrodes |9l [2H or directly through the Coulomb interac- 
tion iTTQl |T7l [22l [23ll . 

The Coulomb interaction potential for two trapped ions, a 
and b, with charges q a and q b in potential wells separated by 
a distance so is given by 
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Here x a and x b are the displacements of the ions from the ex- 
ternal potential minima and £o is the permittivity of free space. 
The first term is constant and does not affect the dynamics. 
The second term represents a steady force between the ions 
that displaces them slightly; if necessary, it can be counter- 
acted with additional potentials applied to nearby electrodes. 
The terms proportional to x\ and x\ represent static changes in 
the trap frequencies that could also be compensated with po- 
tentials applied to nearby electrodes. The term proportional to 
x a x b represents the lowest-order coupling between the ions' 
motions. For small deviations, x' a and x f b , from equilibrium, 
the coupling is 
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and a, a\ b, and represent the harmonic oscillator lowering 
(a, b) and raising (a\ tf) operators, mi and c% are respec- 
tively the ion masses and motional frequencies, fi is Planck's 
constant divided by 271, and we have neglected fast-rotating 
terms. Minimizing the distance, so, between ions is crucial, 
because for fixed coot, the coupling rate scales as £l QX 1 Ao- 
When GX) a = ti>0b — (the resonance condition), we find 

a\t) = exp(/0)oO(^ t (0)cos(^exO-^ t ( ) sm (^exO)( 3 ) 
b\t) = exp(/o)o0(^ t (0) cos (^ex0 -/a t (0)sin(flt ex r)). 

At time t = T ex = ;r/2£l ex , the operators have changed roles up 
to a phase factor, and the oscillators have completely swapped 
their energies, regardless of their initial states. At t = 2r ex , the 
energies have returned to their initial values in each ion. The 
mean occupation, (eft a), of ion a as a function of time exhibits 
oscillations with period 2T ex . 

Figure [T] shows a micrograph of our surface-electrode 
trap 1241 , constructed of gold electrodes, 8 jUm thick with 5- 
fim gaps, electroplated onto a crystalline quartz substrate. The 
trap can produce two potential minima at a height do = 40 fim 
above the surface and separated by so = 40 jim along the x 
axis. Each potential well confines a single 9 Be + ion with an 
axial (parallel to x) motional frequency of COo/ (2k) ~ 4 MHz 
and a barrier between wells of ~ 3 meV. Pseudopotential con- 
finement in radial directions (normal to x) is accomplished 
with a peak potential of ~ 100 V at 170 MHz applied to 
the radio-frequency electrodes, yielding radial frequencies of 
~ 22 MHz. By applying static potentials to the d.c. elec- 
trodes, we can independently vary the separation between 
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FIG. 1: Micrograph of the ion trap, showing radio-frequency (RF) 
and d.c. electrodes, and gaps between electrodes (darker areas). The 
lower part of the figure indicates the simulated potential along the 
trap x axis. Two trapping wells are separated by 40 jim, with ion 
positions marked by red spheres. The d.c. electrodes are sufficient 
to control the axial frequency and the position of each ion indepen- 
dently. Here both frequencies are ~ 4 MHz and the potential barrier 
between the two ions is ~ 3 meV. 

the ions and the curvatures of the two trapping wells. In 
this way, the ion axial motional frequencies can be brought 
into or out of resonance, allowing a tunable interaction. For 
coo /2n = 4.04 MHz, we predict that T ex = 162 /is, where we 
have included a 2% correction in £2 ex owing to the metallic 
electrodes beneath the ions (Methods Summary). 

With currently achieved size scales in ion traps, the 
Coulomb interaction is relatively weak, so low ion heating 
rates and stable trapping potentials are essential. Heating 
rates can be suppressed by operating at cryogenic temper- 
atures (25j [26), so that the direct Coulomb coupling rate 
can exceed the heating rate. In this work, the trap elec- 
trodes and surrounding vacuum enclosure are cooled to 4.2 K 
with a liquid helium bath cryostat. With similar versions 
of this trap at (Oq/2tz = 2.3 MHz, heating rates expressed 
as d(n)/dt (where n denotes the quantum number of mo- 
tional Fock state \n) and (n) is its expectation value) were 
observed to be as low as 70 quanta per second, consistent 
with the results of ref . |26| However, for the experiments de- 
scribed here (coo/271 w 4.0 — 5.6 MHz) the heating rate was 
~ 500 — 2, 000 quanta per second and varied between the two 
wells. We observed d(n)/dt 1/cOq in this trap, in agree- 
ment with previous reports | 25 -27], so large values of C0o are 
beneficial. The use of 9 Be + , the lightest of the commonly 
trapped atomic ions, is an advantage here, because for given 
d.c. trapping potentials the heating rate should remain un- 



changed while £2 ex oc m -1 / 2 . Cryogenic operation decreases 
the background gas pressure to negligible levels, such that ion 
loss rates due to collisions with background gas are smaller 
than one per day. 

A signature of coupling between the ions is the splitting 
between the two axial normal mode frequencies. As the trap 
potential is tuned into the resonance condition, this splitting, 
5/, reaches a theoretical minimum 5/ = £l ex / % = 3.1 kHz. A 
plot of the mode frequencies will therefore show an avoided 
crossing. We measure the mode frequencies by applying a 
nearly resonant oscillating potential pulse to one of the trap 
electrodes. We then illuminate both ions with laser radiation 
resonant with the 2 Si/ 2 — 2 ^3/2 cycling transition at 313 nm. 
A decrease in the resulting fluorescence indicates that a mode 
of the ions' motion has been resonantly excited. For pulse 
lengths T p ^> 1/5/, we resolve the two modes (Fig. [2^). We 
sweep the trapping wells through resonance by varying the 
static potentials that are applied to the trap electrodes. A plot 
of the resulting mode frequencies, determined as above, is 
given in Fig.[2j), c, showing a minimum of 5/ = 3.0(5) kHz, 
in agreement with theory. 

To demonstrate coupling at the level of a few motional 
quanta, we observe the exchange of energy between the two 
ions as follows. The ion motional frequencies are initially 
detuned by 100 kHz, which is much greater than £2 ex /(27r), 
effectively decoupling the ions' motions. The ions are then 
simultaneously illuminated with a laser detuned by —10 MHz 
from the 2 Si / 2 — 2 P?> /2 cycling transition, cooling them into a 
thermal state at the Doppler limit with mean occupation (n) = 
2.3(1). Subsequently, ion a is cooled to (n a ) =0.35(2) by sev- 
eral cycles of stimulated Raman cooling with the 1 1} = \F = 
2,m F = -2) and 1 1) = |1, -1) hyperfine states (2H). The Ra- 
man beams are counter-propagating and oriented at 45° rel- 
ative to the x axis. At this point the potentials are brought 
into resonance ((Dq/2ti = 4.04 MHz) within an interval (9 /is) 
short in comparison with T ex but long in comparison with the 
axial oscillation period. They remain on resonance for a time 
T, allowing energy to transfer between the ions. After a time T, 
the potentials are adiabatically returned to their off-resonance 
values and we determine the mean quantum number, (n a ), in 
ion a by observing asymmetry between the red and blue mo- 
tional sidebands of the | |)-to-| t) hyperfine Raman transi- 
tion l28l 

As seen in Fig. [3] energy exchanges between the ions dur- 
ing an interval T ex = 155(1) /is. The 5% disagreement be- 
tween the measured and predicted (162 /is) values for T ex is 
probably due to uncertainty in the ion separation, so (even a 
l-/im uncertainty would account for the disagreement). The 
first maximum of (n a ) corresponds to the cooling limit of 
ion b = 0)) =2.3(1) quanta). The underlying linear 

growth in (n a ) corresponds to heating of the ions at a rate of 
h = 1, 885(10) quanta per second (Methods Summary). 

As a final experiment, we demonstrate energy exchange at 
approximately the single-quantum level. Ideally, the experi- 
ment takes the following form. The ions are tuned to the res- 
onance condition throughout and are initially Doppler cooled. 
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FIG. 2: Motional spectroscopy of two coupled ions near the avoided 
crossing, a, Decreases in collected fluorescence occur at values of 
excitation frequency, / r f, corresponding to the ion mode frequencies. 
With t p = 960 /is, the splitting on resonance is resolved, b, c, Mode 
frequencies (b) and mode frequency splitting, Sf (c), for the axial 
normal modes of two ions separated by 40 jUm. Error bars are smaller 
than the size of the points. The data were acquired over a 1-h period, 
and slow variations in ambient potentials gave rise to the fluctuations. 
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FIG. 3: Energy swapping between two ions in separate trapping po- 
tentials at the level of a few quanta. The mean occupation, (n a ), of 
ion a is plotted with error bars (s.e.m.) for various durations, T, that 
the ion motional frequencies remain on resonance. The blue curve 
represents a fit to theory with four free parameters: the two initial 
mean quantum numbers, the exchange time and the heating rate. En- 
ergy exchanges between the ions at 155(l)- J us intervals. The linearly 
increasing trend in (n a ) is due to ion heating at a rate of 1,885(10) 
quanta per second. Uncertainties represent standard errors of the fit 
parameters. 
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FIG. 4: Motional exchange between two ions in separate trapping 
potentials at approximately the single-quantum level. The probabil- 
ity, P(\ t)a)» of measuring ion a in spin state | t)« at the en d of the 
experimental sequence is plotted with error bars (s.e.m.) against the 
time, T, for which the ions interact. P(| t)<a) oscillates with period 
2r ex = 437(4) ^s as a quantum exchanges between the ions. 



Ion a is Raman-cooled, sympathetically cooling ion b and 
thereby preparing the state |0) fl | l) a \0)b' To create a single 
motional quantum, we drive ion a with a blue- sideband Ra- 
man n pulse (of duration 10 /is, which is much less than T ex ), 
creating the state |l) fl | t)a|0)&. The system oscillates between 
1)a\0)b and |0) fl | t)a\l)b with period 2 T ex . After a time T, 
we drive ion a with another blue-sideband n pulse, condition- 
ally flipping the spin from |l) fl | f) a to |0) fl ||) a , dependent on 
the presence of a motional quantum in ion a. The final internal 
state probability will be given by P(| X)a) (t) = sin 2 (£2 ex T). In 
practice, contrast in the oscillations (Fig. [4]) is significantly re- 
duced by incomplete cooling, motional decoherence and de- 
coherence due to imperfect Raman sideband pulses. We es- 
timate that ion a is cooled initially to (n a ) = 0.3(1). Al- 



though we were unable to measure the initial temperature of 
ion b directly, comparison of the contrast and temporal behav- 
ior of our exchange data (Fig. [4]) with simulations indicates 
that (rib) i$ 0-6. Motional decoherence results from heating 
and from trap frequency instability over the time required to 
acquire the data. Raman sideband pulses suffer from varia- 
tions in laser intensity and fluctuations in the sideband cou- 
pling caused by thermal spread in the y and z motional states 
(Debye-Waller factors ifTOl ). For C0o/2n = 5.56 MHz, we ob- 
serve oscillations with period 2r ex = 437(4) /is. The 2% dis- 
agreement with the prediction (447 /is) of equations [2] and [3] 
is probably due to uncertainty in the ion separation and the 
difficulty of maintaining the exact resonance condition. 
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Significant improvements in coupling fidelity seem to be in 
reach with current technology. It should be possible to im- 
prove Raman laser intensity stability, and Debye-Waller fac- 
tors from the y and z motion can be eliminated by proper 
choice of beam directions l28l . We believe that our motional 
frequency instability can be alleviated by ensuring that the 
trap surface is free of charged contaminants. Faster exchange 
can be achieved by scaling down trap dimensions, but this puts 
a premium on reduced motional heating. Although the work 
presented here uses the axial mode to couple the ions, it may 
prove advantageous to use the radial modes in future experi- 
ments because of the lower heating rates associated with their 
higher frequencies. As an example, radial mode frequencies 
of ~ 30 MHz are routinely achieved in the apparatus, com- 
pared with axial frequencies < 10 MHz. 

These results could lead to several possible applications in 
quantum state engineering and spectroscopy. For example, 
from the motional state |l) a |0)& the state of the ions at time 
t = T ex /2 is the Bell state (\l) a \0) b + i\0) a \l) b )/\/2. Trans- 
ferring this state onto the ions' internal states with sideband 
pulses would create an entangled spin state, even between ions 
of dissimilar species. This could be used as an entangled- 
pair factory in the scheme of refs 1 101 181 with the advantage 
over previous schemes Q that the ions are already in separate 
wells, ready for distribution to separate locations. The cou- 
pling could be used to read out the state of one ion species 
with another, an ability useful for error correction protocols 
and for quantum logic spectroscopy fT9l . When the harmonic 
wells are not in resonance, the spin state of one ion can be 
read out without destroying the state of the other, so schemes 
for weak or quantum non-demolition measurements that use 
either Kerr-type nonlinearities or quantum logic may become 
feasible Q. Hybrid quantum systems, incorporating similar 
interactions to couple a trapped ion to other quantum devices, 
could serve as a means of transferring quantum information 
between different qubit implementations in a future quantum 
network. For example, a trapped ion could act as a quan- 
tum transformer between a superconducting qubit 031 and 
a photonic qubit (221 EI • The sympathetic cooling through 
exchange might also be used to cool neutral molecules l3T1l . 

METHODS SUMMARY 

The shielding factor, /3, represents the ratio of the exchange 
rates £l QX with and without the presence of the trapping elec- 
trodes. To a good approximation, we ensure that at the 
motional frequencies all trap electrodes are held at ground. 
Therefore, assuming that gaps between the electrodes are neg- 
ligible, the shielding factor can be calculated with the method 
of images. The result is 

1 / 3(soM)) 5 (s /d ) 3 \ 

2 V (4 + (soM)) 2 ) 5/2 (4 + (*oM)) 2 ) 3/2 J 

which reaches a maximum of j3 = 1.018 at so = do. 

The evolution of (n a ) under the influence of equation [T] 
including heating effects and assuming that both ions begin in 



a thermal state with mean quantum numbers n a o and n b o, can 
be predicted with a Langevin equation 0. On resonance, the 
evolution is 

(n a ) = n a0 cos 2 (£l ex t) + n b0 sin 2 (£l t xt) +nt (5) 

where h represents the mean of d(n a )/dt and d(n b )/dt for 
uncorrected noise sources. 
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